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Abstract. A quantum mechanical description of particle propagation on the discrete 
spacetime of a causal set is presented. The model involves a discrete path integral 
in which trajectories within the causal set are summed over to obtain a particle 
propagator. The sum-over-trajectories is achieved by a matrix geometric series. For 
causal sets generated by sprinkling points into 1+1 and 3+1 dimensional Minkowski 
spacetime the propagator calculated on the causal set is shown to agree, in a suitable 
sense, with the causal retarded propagator for the Klein-Gordon equation. The 
particle propagator described here is a step towards quantum field theory on causal 
set spacetime. 

PACS numbers: 03.65.Pm 04.60.Nc 
1. Introduction 

Causal set theory is an approach to quantum gravity in which spacetime is fundamentally 
discrete. Spacetime events are represented by elements of a causal set — a locally finite, 
partially ordered set in which the partial order represents the causal relationships 
between events in spacetime. The reader is directed to [U [2] for detailed introductions, 
motivations and further references. 

Modelling spacetime as a causal set is a radical departure from the usual model of 
spacetime as a Lorentzian manifold. As such we expect the concepts used to describe 
physics on a Lorentzian manifold will not be directly applicable to describing physics 
on a causal set. In particular the differential equations of motion used in the continuum 
are clearly not appropriate (at least in an unmodified form) if spacetime is modelled as 
a discrete structure. The discrete path integrals presented here (see Section [3]) are an 
example of applying continuum ideas to a causal set in a natural way. 

The appropriate way to model matter on a causal set is an important physical 
question which, if answered, could provide testable predictions from the theory. The best 
description of matter is currently quantum field theory and we expect that each reader 
will have his or her own mental picture for interpreting quantum field theory — a picture 
based on fields, waves, particles or some combination of these. In this paper we take 
the viewpoint that matter is modelled as point particles (rather than fields) propagating 
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according to quantum mechanical laws. This is in keeping with the propagator approach 
to quantum field theory pioneered by Feynman (See [10J). A particle description rather 
than a field description is also, perhaps, more in keeping with the "discrete" causal set 
philosophy. 

In Section [3] we perform quantum mechanical path integrals (or rather path sums) 
on a causal set. Due to the discrete nature of a causal set these can be defined as 
actual sums over all particle trajectories (which move forwards in time) leading from 
one spacetime event to another. A particle propagator is defined by summing amplitudes 
assigned to particle trajectories which move forwards in time. We show that, for causal 
sets generated by sprinklings into 1+1 and 3+1 dimensional Minkowski spacetime, this 
propagator gives (in a suitable sense) the Klein-Gordon retarded causal propagator from 
quantum field theory. This describes the propagation of free scalar particles provided 
there is no actual pair production in the vacuum [TT]. 

Applying the path integral formalism to discrete spacetime has been considered 
before. Jacobson [3] and Gudder [5] have separately considered path integrals on a 
hypercubic lattice and similar path integral ideas for causal set spacetime have been 
considered by Meyer jH] and Foster and Jacobson [I]. Previous approaches to describe 
matter on discrete spacetime have often relied on finite-differencing the continuum 
equations (see, for example, [7j). If we are developing a fundamentally discrete theory, 
however, the continuum description should emerge from the underlying theory. The 
path integrals presented here are an example of an underlying discrete theory leading 
to a continuum propagator. 

2. Definitions 

A causal set (or causet) is a locally finite partially ordered set. This means it is a pair 
(C, :<) with a set C and a partial order relation ^ defined on C that is 

• Reflexive: x ^ x, 

• Antisymmetric: x <y <x implies x = y, 

• Transitive: x <y < z implies x -< z, 

• Locally finite: \{z G C\x -< z -< y}\ < oo for all x, y G C, 

for all x,y,z G C . Here \A\ denotes the cardinality of a set A. We write x -< y to mean 
x <y and x ^ y. 

The set C represents the set of spacetime events and the partial order ^ represents 
the causal order between pairs of events. If x -< y we say u x precedes y" or "x is to the 
causal past of y" . 

A chain in a causal set C is a subset in which any two elements are related. If the 
chain is finite it is a sequence of distinct elements Vi , . . . , Vi n each of which precedes the 
next element in the sequence: Vi -< v ^ -< v j 2 -< . . . -< v,- ln . The length of the chain is n. 

A link (or covering relation) between u, v G C is a relation u -< v such that there 
exists noweC with u -< w -< v. We say u and v are nearest neighbours (or v covers 
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u) and write u -<* v. 

A path is a subset P C C which is a maximal (or saturated) chain. This means P 
is a chain with no element w G C — P such that u -< w -< v for some u,v £ P. If the 
path is finite it is a sequence of distinct elements Vi , . . . , Vi n each of which is linked to 
the next element in the sequence: v io -<*v il ^*v i2 -<*... -<*v in . The length of the 
path is n, the number of links. 

2.1. Adjacency matrices 

A useful way to represent a finite causal set is by its adjacency matrix. We now define 
two types of adjacency matrices that will be used to compute the path integrals in 
Section |3j We begin by looking at finite causal sets as this allows us to work with finite 
square matrices. In Section [4] we extend the work to general causal sets. 

Before defining the adjacency matrices we must label the elements in our causal 
set. For a finite causal set C with p elements we choose a labelling v%, . . . ,v p for the 
elements of C . There are then two adjacency matrices of interest: Ac and Ar. These 
are p x p matrices defined by 

(*ch ■= | I 
(Afi)ij := | J 

These are both zero on the main diagonal and, from the definition of a causal set, 
the labelling can always be chosen to ensure they are both strictly upper triangular 
matrices (in which case the labelling is called a natural labelling). Within graph theory, 
the matrix Ac is called the transitively closed adjacency matrix and Ar is called the 
transitively reduced adjacency matrix. Sorkin [17] has termed Ac the causal matrix and 
Ar the link matrix. 

Powers of these matrices have the following useful properties [Hj: 

(Ac)ij = The number of chains of length n from V{ to Vj, (2-3) 
(^S)ij = The number of paths of length n from V{ to Vj. (2-4) 

2.2. Sprinklings 

There exist causal sets which bear little resemblance to the smooth Lorentzian manifolds 
of general relativity. It is hoped that in a full theory of causal set dynamics the causal 
sets which are "manifold-like" on large scales will emerge dynamically from the theory. 
For the time being, however, if we are interested in comparing the results of calculations 
performed on a causal set with calculations on a Lorentzian manifold we must ensure 
by hand that the causal set is like the manifold. 

We do this by performing a sprinkling (see [9] for more details) into a (^-dimensional 
Lorentzian manifold (M, g). To perform a sprinkling points are placed at random within 



if Vi -< vj 
otherwise, 

if Vi ^* Vj 
otherwise. 
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M using a Poisson process with sprinkling density p. This means that the points are 
randomly distributed throughout M with n points being sprinkled into a region of d- 
volume V with probability 

(pV) n 

Prob(n points in region of d-volume V) = ^-^-e~ pv . (2.5) 

n! 

This ensures the expected number of points in a region of d- volume V is equal to pV. The 
sprinkling generates a causal set whose elements are the sprinkled points. The causal 
set's partial order relation is "read off" from the manifold's causal relation restricted 
to the sprinkled points. If the Lorentzian manifold satisfies a strong enough causality 
condition (e.g. if it's strongly causal) then we expect the manifold to well-approximate 
the sprinkled causal set. 



We note that setting n = in (2.5) gives the probability that two points x -< y in 
M have no sprinkled points causally between them as e~ pv( - x ~ y \ Here V(x — y) is the 
d- volume of the causal interval in M between x and y. 

3. Path integrals on causal sets 

A classical model for particle propagation is the swerves model described in [H] , As 
acknowledged in that paper, however, a fundamental description should incorporate the 
principles of quantum mechanics. 

In the continuum the propagation of particles is described by a propagator K(x, y) 
which can be interpreted as the quantum mechanical amplitude (or just amplitude) for 
a particle to travel between two spacetime points x and y (see, for example, [TOj chapter 
6]). The path integral approach to quantum mechanics [T21 H3] allows the propagator 
to be calculated by assigning amplitudes to each possible particle trajectory from x 
to y. Summing these amplitudes over all possible trajectories from x to y gives the 
propagator K(x,y). As we shall show, in the discrete spacetime of a causal set this sum 
over trajectories can be simply performed. 

To define path integrals on a causal set we have to make two choices: which 
trajectories to sum over and what amplitudes to assign to each trajectory. In this 
paper we'll consider two types of trajectories within the causal set: either all possible 
chains between two elements or all possible paths between two elements. The amplitude 
assigned to each trajectory is then defined in terms of two basic amplitudes: a and b. 

The constant a is the amplitude for the particle to 'hop' once along the trajectory 
from one element to the next. The constant b is the amplitude for the particle to 'stop' 
at one element of the trajectory (the initial and final elements are not regarded as stops). 
The amplitude for the whole trajectory is then the product of the amplitudes for each 
hop and each stop. For a chain or path of length n (so there are n hops and n — 1 
intermediate stops) this compound amplitude is therefore a n b n ~ x . 

In Section |3.2| we identify values of a and b which, for causal sets generated by 
sprinkling into 1+1 and 3+1 dimensional Minkowski spacetime, give (in a sense to be 



defined) the causal retarded Klein-Gordon propagator. These values (equations (3.31) 
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and (3.44)) depend on the mass of the particle and on the spacetime volume assigned 
to each causal set point. 

As in Section 2A_ we begin by working with finite causal sets. How the ideas extend 
to arbitrary causal sets is presented in Section |4j For a causal set with p elements we 
define a p x p matrix $. 

If we sum over chains we define 

$ := aA c . (3.1) 

If we sum over paths we define 

$ := aA R . (3.2) 

The total amplitude to go from Vi to Vj along a trajectory of any length is then Kij 
where K is the p x p matrix 

oo 

K :=/ + $ + b<t> 2 + 6 2 $ 3 + . . . = I + b 71 - 1 ^. (3.3) 

rc=l 

Here matrix multiplication is used and I is the p x p identity matrix. We have included 
the / matrix to ensure = 1 for alH = 1, . . . ,p although this is just convention. 

Each term in this sum is the contribution to the total amplitude from chains or 
paths of a particular length. As an example, 

v 

H$% = $>$Ai, (3-4) 

k=l 

is the sum (over all intermediate positions Vk) of the amplitudes for all length two 
trajectories Vi to Vk to Vj. 

Since the causal set is finite and the trajectories move only forwards in time the 
sum in (3.3) terminates and we have 

K = / + $(/- (3.5) 

This matrix inverse is simple to perform if the causal set has been labelled to ensure 
$ is strictly upper triangular. In this case the upper triangular matrix 7 — 6$ can be 
inverted using elementary row operations. 

The question we now face is whether there are values of a and b such that K^, 
for a causal set generated by a sprinkling into Minkowski spacetime, is approximately 
equal to the continuum propagator in Minkowski spacetime. We address this in the 
next section. 



3.1. Expected values 

Causal sets generated by sprinklings into a Lorentzian manifold will not be identical 
because the particular points that are sprinkled are chosen randomly. The value of the 
propagator calculated from one sprinkled causal set therefore depends in detail on the 
particular causal set that is generated. This randomness in the sprinkling process means 
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we are interested in the expected value (and variance) for the propagator calculated by 
averaging over all possible sprinkled causal sets. 

We shall look at sprinklings into (^-dimensional Minkowski spacetime. Firstly we 
shall fix two points x and y in Minkowski spacetime. We then sprinkle a causal set. 
There is zero probability that the sprinkled causal set will contain x and y so we then 
add x and y to the sprinkled causal set. In this manner, averaging over all sprinkled 
causal sets (with a fixed sprinkling density p), we shall calculate the expected number 
of chains and paths from x to y and the expected value of the propagator between x 
and y. 



3.1.1. Summing over chains To calculate the expected number of chains between two 
points x and y in ci-dimensional Minkowski spacetime we define 

v{x-y):=[ l \ X ~ y (3.6) 
v yj 1 otherwise. v ; 

Translation invariance of Minkowski spacetime ensures v is only a function of the 
separation x — y. The expected number of chains of length one from x to y is given by 

C\(x-y) = v{x-y). (3.7) 

The expected number of chains x -< z\ -< . . . -< z n -\ -< y of length n > 1 is given by [6] 

C n {x - y) := p 11 - 1 J - J d d Zl ■ ■ ■ d d z n ^v{x - z x )v(z x -z 2 )--- v{z n - X - y). (3.8) 

The expected value for a propagator which sums amplitudes assigned to chains is thus 

oo 

K c (x -y):=J2 a n t> n ~ l C n {x - y). (3.9) 

n=l 

This satisfies the integral equation 

— y) = av[x — y) + abp J d d z v{x — z)Kq{z — y). (3.10) 

3.1.2. Summing over paths To calculate the expected number of paths between two 
sprinkled points x and y in <i-dimensional Minkowski spacetime we define 

^-^)-| otherwise . (3- 11 ) 

where V{x — y) is the <i-dimensional Minkowski spacetime volume of the causal interval 
between x and y. The expected number of paths of length one from x to y is given by 

Pi(x - y) = p P (x - y). (3.12) 

The expected number of paths x -<* z\ -<* ... -<* z n -\ -<*y of length n > 1 is given 
by[T6] 



P n {x -y) = p 



n-1 



J " J ' ' ' & d Zn-iV<p{x ~ z x )p p (z x - z 2 ) ■ ■ ■ p P {z n -i - y). (3.13) 
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The expected value for the propagator which sums over paths is thus 

oo 

K P (x-y) :=J2 anbn ~ lp n(x-y)- (3-14) 

71=1 

This satisfies the integral equation 

(, 15 , 

3.2. Minkowski spacetime propagators 

In (i-dimensional Minkowski spacetime the causal relations are particularly simple. For 

x = (x ,x), y = (y ,y) we have: 

x r< y x < y and (x - y ) 2 > (x - y) 2 . (3.16) 

The scalar propagators in (i-dimensional Minkowski spacetime quantum field theory 
are Green's functions of the Klein-Gordon equation: 

{U + m 2 )K^\x-y) = 5 d {x-y). (3.17) 

Here m is the mass of the particle, 5 d is the (i-dimensional Dirac delta function and we 
choose units with h = c = 1. The d'Alembertian is given by 

Jf_ JP_ JP_ _ d 2 
dx 2 dxi 2 dx 2 2 "' dx d ^i 2 ' 

To obtain the propagator Km (% ~ y) explicitly we define the Fourier transform by 

f(p) := J d d xf(x)e^, f(x) = j^y d j d d pf(p)e-^, (3.19) 



where px := po^o — p ■ x. Using this we can solve equation (3.17) to obtain 



K£>(p) := -^-^ ~ 2 , (3.20) 

Po ~ P m 

so 

-ip{x-y) 



K%\x-y) :=-t^ / d W^- , . (3.21) 



e 

(2ir) d J ^ Po — p 2 — m 2 

( d )] _ T\/Td~2 



The mass-dimension of this propagator is [K m ] = M a 

The integrand in these expressions contains poles so the propagator in Minkowski 
spacetime is only uniquely defined if we specify a contour of integration (or equivalently 



boundary conditions for the solution of (3.17)). The retarded causal propagator is the 



Green's function obtained by avoiding the poles at po = ^^/p 2 + m 2 by two small 
semi-circles in the upper-half po complex plane. This is equivalent to 

1 f p -ip(x-y) 

K^ix-y) := lim -— — - / d d p- -. (3.22) 

mK yJ e^o+ (2vr) d J y (p + it) 2 - p 2 - m 2 V ; 

This can be calculated explicitly in different dimensions. We consider the expressions 
in d — 2 and d — 4 dimensions. 
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In 1+1 dimensional Minkowski spacetime we have: 

K V>(x-y) = [ ^ o(mr - } \ X ^ y (3.23) 

m K y! y otherwise, v 7 

^i 2) (p) = -7 , .J 2 a- (3-24) 

(Po + «0 2 -Pi-m 2 

Here r xy = ■>/ (x — ?/o) 2 — (#i — 2/i) 2 is the proper time from x to y and J is the zeroth- 
order Bessel function of the first kind. 

In 3+1 dimensional Minkowski spacetime [15J: 

f 1 X( T 2 \ m Ji(mr xy ) -r 

K^\x-y) = \ ^ [xy) 4 - ~ y (3.25) 
I otherwise, 

K ( *\p) = -7-—^^ 2" ( 3 - 26 ) 

(Po + ze)^ — pr — m 



Here r xy = (xq — yo) 2 — (x — y) 2 is the proper time from x toy and J\ is the first-order 
Bessel function of the first kind. 

We keep the e terms in the Fourier transformed propagators to ensure the poles 
are avoided in the correct way. It is understood that the e — > + limit is taken at the 
end of the calculations. It is these position-space propagators we wish to reproduce on 
a causal set. 



3.2.1. 1+1 dimensional Minkowski spacetime It turns out that in 1+1 dimensions the 
propagator requires us to sum over chains. Fourier transforming (3.10) the integral, 
being a convolution, becomes a product and we have 

K c (p) = au(p) + abpu(p)K c (p), (3.27) 

or 

av(p) 



1 — abpu(p) 

In 1 + 1 dimensions the function v{x — y) has Fourier transform 

2 



u(p) 



(2) 

since v = 2K ( q> (and using ( |3.23[ ) and ( |3.24[ ) with m — 0). Substituting into ( J3.28[ ) we 
have 

In 

2a 



(po + ie) 2 - p\ 



2 ' 



(3.28) 



(3.29) 



K G {p) 



(p +ie) 2 -pl 



1 + 



2abp 



{po+ie) 2 -pf 



(p + it) 2 - p\ + 2abp 



Equating this to equation (3.24) we find 

1 m 2 
a — -, b = — - 



(3.30) 



(3.31) 



2' P 
are the correct amplitudes. 

In 1+1 dimensions \p\ = M 2 so the amplitudes assigned to the chains are 
dimensionless: [a] = [b] = [a n 6 n_1 ] = 1 (for n = 1, 2, . . .). 
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3.2.2. 3 + 1 dimensional Minkowski spacetime It turns out that in 3+1 dimensions 
the propagator requires us to sum over paths. Fourier transforming (3.15) we get 

ap p {p) 



(3.32) 



abpp p (p) ' 

In 3+1 dimensional Minkowski spacetime V(x — y) := ^r^ y is the volume of the causal 
interval between x and y. We therefore have 

e -pV{x-y) 



M x - y) 



o 



e 24P T xy if x ^ y 

otherwise. 



The function 



if z < 



(with a real constant c > 0) satisfies 

lim f p {z,c) = ^—rdiz), 
p^oo £y c 

where 5(z) is the Dirac delta function. We therefore see 



lim Jp~ n P [x — y) = lim 

p^oo p^oo 



otherwise 



xyJ 







if # d: y 

otherwise. 



Setting m = in (3.25) and (3.26) we see that in 3+1 dimensions we have 
K^\x-y) = 







2 ) 

xyJ 



if x ^ y 



otherwise, 



^ 4) (P) 



1 



(p + ^e) 2 - P 



These imply, taking Fourier transforms and using (3.36), that 

\/P Pp(P) = 27T- 

p^oo 

Setting 

B 



lim y/p fi p {p) = 27i^K^\p) = -27rv / 6 7 — -y 2 

2 (po + «e) -p 



a = Ay/p, 



P 



(3.33) 
(3.34) 
(3.35) 
(3.36) 

(3.37) 
(3.38) 

(3.39) 
(3.40) 



where A and B are (possibly dimensionful) constants independent of p we substitute 



into (3.32) to get 

Kp(p) 



l-AB^pp p (p) 
As p tends to infinity this becomes 

AC 



lim Kp(p) 

p^oo 



AC 



ABC 



{po + ief -p> + ABC' 1 



(3.41) 



(3.42) 



where C := 27rv6. 
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Equating this with (3.26) we have 

AC = 1, B = -m 2 , (3.43) 

so 

^ m * (3.44) 



27^' p ' 

are the correct amplitudes. 

In 3+1 dimensions [p] = M 4 so the amplitudes assigned to the paths have mass- 
dimension M 2 : [a] = M 2 , [b] = M~ 2 , [a 71 ^' 1 } = M 2 (for n = 1,2, . . .). 

3.3. Fluctuations 



The basic amplitudes, a and b, given in (3.31) and (3.44) ensure that the expected 



value of the propagator, averaged over causal sets generated by sprinkling into 1+1 and 
3+1 dimensional Minkowski spacetime, equals the causal retarded propagators for the 
Klein-Gordon equation (with the infinite density limit being taken in the 3+1 case). 
The propagator values calculated on any particular sprinkled causal set will, however, 
fluctuate randomly away from the expected values. An important question is how 
the size of these fluctuations vary as the density of the sprinkling increases. If the 
fluctuations increase as the density increases (i.e. in the continuum limit) then the 
model clearly needs to be changed. An example of such a situation is given in [T7] 
where a causal set approximation to the d'Alembertian has growing fluctuations in the 
large density limit. This is resolved by modifying the d'Alembertian approximation by 
performing a sum over points contained within a new non-locality scale of the original 
point. 

One way to investigate the question of how the fluctuations depend on the sprinkling 
density is through numerical simulations. These involve generating causal sets by 
sprinkling points into an interval of Minkowski spacetime. We then calculate the 
propagator on the generated causal set and examine the fluctuations for different density 
sprinklings. The fluctuations depend on the sprinkling density p as well as on the size 
of the particle's mass m relative to p. The simulations discussed here use causal sets 
generated by David Rideout's Cactus computer code (briefly described in [IS]). 

For sprinklings into 1+1 dimensional Minkowski spacetime and for m 2 <C p 
preliminary results suggest the fluctuations decrease as the density increases. 

For sprinklings into 3+1 dimensional Minkowski spacetime the propagator's 
expected value equals the continuum propagator only in the infinite density limit. A 
realistic sprinkling density would be larg^f] but finite. Investigating the behaviour of 
the propagator for large finite sprinkling densities is difficult to do through numerical 
simulations. This is because current simulations cannot produce enough sprinkled points 
to ensure a large density over a large spacetime volume. 

| Equivalently we can use a small density but sprinkle into a Minkowski spacetime interval of large 
volume. 
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High density sprinklings can be achieved by sprinkling a moderate number of points 
into a small volume. In this case, however, the behaviour of the propagator is only 
investigated within such small volumes. In 3+1 dimensional Minkowski spacetime 
the preliminary results from such simulations suggest, for m 2 <C ^Tp, the expected 
value of the propagator for large densities well-approximates the continuum propagator. 
In addition we may very preliminarily say the fluctuations away from the continuum 
propagator are small for large density. These conclusions hold only for the small volumes 
considered in the simulations and further work is needed to check whether this behaviour 
holds for large density sprinklings into large spacetime volumes. Of course determining 
the rate at which the fluctuations decrease remains to be done. 

We note here that the constraints on the mass of the particle are easily satisfied 
for realistic sprinkling densities. Assuming a Planckian sprinkling density in the 3+1 
case we let p be the inverse of the Planck 4-volume: p = c 1 j[GK) 2 . The heaviest known 
elementary particle is the top quark with a mass m = 174.2 ± 3.3GeV |19j. In Planck 
units, so p = 1, this is m « 1.4 x 10~ 17 . This certainly satisfies m 2 / \Jp~ ^2x 1CT 34 <C 1. 
If the masses of particles modelled on the causal set are of order the masses of known 
elementary particles then the causal set propagator (for a Planckian density sprinkling) 
should well-approximate the continuum propagator. 

4. General causal sets 

We now consider general causal sets — including non-finite causal sets as well as those 
not generated by sprinkling into Minkowski spacetime. 

4-1. Incidence algebra 

An interval (or Alexandroff neighbourhood) in a causal set C is defined to be 

[u,v] := {w E C\ u ■< w ^ v}, (4.1) 
provided u v. 

Since a causal set is locally finite the propagator K(u, v) for any two points u -< v 
can be calculated by applying the methods of Section[3]to the finite interval [it, v). There 
is another way, however, to view the path integral framework which uses the incidence 
algebra of a causal set. 

For a causal set C (not necessarily finite) we denote the set of all intervals by 

Int(C) := {[u,v]\u,v G C, u ■< v}. (4.2) 

The incidence algebra [8j of C over C, denoted I(C,C), is then the associative 
algebra of all functions 

/ : Int(C) -> C, (4.3) 

with multiplication defined by 

f*g(x,y):= ^ f(x,z)g(z,y). (4.4) 
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The sum here is finite (so f*g is defined) because the causal set is locally finite. I(C, C) 
is an associative algebra with two-sided identity 

We note that we could use an algebraic field other than C in defining the algebra. The 
adjacency matrices, which were defined only for finite causal sets, generalize to the 
algebra elements 

, 1 if x -< y ,. „ N 

A ^v)--={ o other wL, (4 ' 6) 

A R {x,y):={ 1 if f ^ V (4.7) 
* otherwise. v ; 

Powers of these functions, under *, satisfy: 

Ag(x, y) = The number of chains of length n from x to y. (4.8) 

A T ft(x, y) = The number of paths of length n from x to y. (4-9) 

The path integral work in Section [3] can be done using the incidence algebra. 
Phrasing the method this way rather than using adjacency matrices allows us to work 
with infinite causal sets without restricting to a finite sub-causal-set. For finite causal 
sets the incidence algebra and adjacency matrix methods are entirely equivalent. 

If we sum over chains we define an element $ of J(C, C) by 

®(x,y):=aA c (x,y). (4.10) 
If we sum over paths we define 

$(x,y):=aA R (x,y). (4.11) 
We then have, in a manner similar to the finite case, that 

K := 5 + $ + &($ * $) + 6 2 ($ *$*$) + ... (4.12) 

is the algebra element for the propagator. Here K(x, y) is the quantum mechanical 
amplitude that the particle will travel from x to y along a trajectory of any length. 
Formally we have 

K = 6 + $ * (5 - b§)-\ (4.13) 

but it is not immediately clear if an inverse of 5 — 6$ exists. Applying Proposition 
3.6.2 in [8] shows, however, that (5 — 6$) _1 exists so K can be written in this form. 
The proposition also shows that K(x, y) depends only on the values of $ for intervals 
contained within [x,y]. 

It is pleasing to note that the choice of modelling spacetime as a causal set enables 
us to define its incidence algebra which can then be used to define path integrals. The 
choice of causal set spacetime fits naturally with the rules governing quantum mechanical 
amplitudes. 
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4-2. Non-sprinkled causal sets 

The basic amplitudes for sprinkled causal sets depend on the sprinkling density p. For 
causal sets not generated by a sprinkling, however, we must make sense of the p that 
appears in a and b. To do this we assume that an arbitrary causal set element is assigned 
a fundamental spacetime volume Vq. When the causal set is generated by a sprinkling 
into Minkowski spacetime with density p we have Vq = 1/p. 

If the causal set leads to a macroscopically 1+1 dimensional spacetime we sum over 
chains and assume [Vo] = M~ 2 . This gives 

a = -, b=-m 2 V . (4.14) 

If the causal set leads to a macroscopically 3+1 dimensional spacetime we sum over 
paths and assume [Vo] = M~ 4 . This gives 

a =^W b= - m2y °- (4i5) 



5. Further work 



One of the motivations for discrete spacetime is to regularize the divergences in quantum 
field theory The propagators presented here could be used to develop quantum field 
theory on a causal set and, in particular, the graph-based language of Feynman diagrams 
appears particularly well-suited for this. This paper deals only with scalar particles and, 
as required for a realistic quantum field theory, we would hope to model spinor and 
vector particles on a causal set — possibly by allowing internal particle states. Another 
interesting question is whether the Feynman propagator can be obtained by a path 
integral on a causal set — possibly by allowing trajectories which can reverse time- 
directioii[§J These remain open questions. 

The perspective used in this paper is that of relativistic quantum mechanics — the 
motion of matter is described upon a fixed spacetime background (the causal set). No 
backreaction has been included in which the structure of the causal set would depend 
on the matter distribution upon it. This is one of the principal difficulties in quantum 
gravity — to obtain a description of matter and spacetime that is able to describe more 
than just quantum matter on a fixed classical spacetime background. This may be 
achieved in the long sought-after "sum over causal sets" version of the theory. 
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§ These paths would be interpreted as particle pair production[10 
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